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In this work we propose a new mechanism for photoproduction inside a neutron star based on 
Parametric Resonance phenomenon as firstly applied to Inflationary Cosmology. Our assumptions 
are based on the pion condensation model by Harrington and Shepard. We show that a huge 
number of photons are created which, on turns, reheats the matter in the star's core. Thus, we 
argue that Parametric Resonance can be effective during a brief period out of an neutron star 
lifetime leading to an anomalous uprising variation of its brightness departing from the black body 
radiation at regularly spaced frequencies. In adition, a time periodic signal is obtained in moderate 
(not exponential) regimes. We argue also that our PR mechanism offers a simple and feasible 
explanation for some recent observations of giant flares from neutron stars. 

PACS numbers: 97.60. Jd,13. 60. L 



(N 
> 

O 

in 
o 



Of 
i 

o 

CO 



X 



I. INTRODUCTION 

It is well known that in Inflationary Cosmology the 
early Universe turns out be depleted of its matter con- 
tent at the end of inflation era. So a mechanism is neces- 
sary in order to fill up the Universe with baryonic matter 
(baryogenesis) , as well dark matter. In current models, 
matter must be generated at a very high temperature, 
after the inflation era, in order to start baryogenesis. A 
mechanism, firstly proposed by Traschen and Branden- 
berger and independently by Dolgov and KirilovaH 
based on Parametric Resonance (PR) phenomenon 
succeed to provide such requirements as a possible way 
to generate matter in pos-inflationary Universe with an 
exponential particle production. Few years later, Kofman 
et al Q showed that this mechanism is effective mainly 
in the so-called preheating phase of the inflation era. 

The simplest model in Field Theory consists of two 
boson fields, namely, a self-interacting source field (f> and 
a resonant (or matter) field x- From the mathematical 
point of view, the source field </> can be obtained as so- 
lution of an ordinary second order nonlinear differential 
equation. Periodic solutions are founded when the field 
(j) is in the oscillatory regime around the minimum of its 
potential (for example, in A</> 4 model). In the simplest 
(spatial homogeneous) case, one looks for periodic solu- 
tions which in turns enter in the ordinary linear (time) 
differential equation of the matter field \ as a periodic 
potential. These linear equations are particular cases of 
the so called Hill's Equation The important point 

for these models is that Hill's Equation presents, in its 
parameter space, instability regions, also named resonant 
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bands. In a resonance band the amplitude of a x solution 
can grow exponentially and, since the particle number is 
proportional to this amplitude one obtains a huge num- 
ber of x quanta 4]. If x is a boson field, the particle 
number occupation with a specific momentum is expo- 
nentially large. On the other hand, Dolgov and Kirilova 
showed that the PR mechanism is not very effective 
(particle production is not exponential) for fermions field 
due to Pauli blocking and the solutions presents only a 
polynomial growth. 

Few years ago, Garcia-Bellido and Kusenko 7] pro- 
posed that the same PR mechanism could be used to ex- 
plain the huge energy coming from "point-like" sources 
in the sky, i.e. the so-called cosmic ray bursts. In order 
to obtain the needed potential, bounded from below, the 
above authors took the collision of two neutron stars (NS) 
which, merging in each other produce a very compact ob- 
ject whose proton matter inside turns out to be supercon- 
ductor and can be modelled by a Ginzburg-Landau po- 
tential ||. In addition, they took electromagnetic field 
as a resonant field interacting with the superconductor 
matter. So, when the system is in the PR regime a huge 
burst of energy is transferred from superconductor mat- 
ter to the electromagnetic field in form of high energy 
photons and these, in principle, could be detected on the 
Earth as cosmic ray bursts. As far as we known this was 
the first application of the PR mechanism in order to 
explain the bursts coming from compact objects in the 
sky. 

Nevertheless we wonder whether PR can be applied, 
not to the spectacular collision of two NS, but to a sin- 
gle star itself in a more moderate energy scale, of course. 
This was motivated, after Walecka's pioneering work || 
on quantum field theory applied to stellar structure. So, 
differently from the above cited researchers, we apply the 
PR mechanism to a single NS, whose part of its matter 
content in the core has supranuclear densities. We are 
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mainly interested in the regime in which matter in the 
star medium is populated with negative pion described 
by a field in a condensate phase ^3] which oscillates 
around the minimum of its potential as required by the 
PR mechanism. This oscillation is a natural consequence 
of the energy excess coming from the formation of the 
pion condensate, that is, the system to turns in a con- 
densate phase must reduces its energy and therefore is 
plausible to suppose that this excess is then transferred 
to the condensate field causing its motion. As in Garci'a- 
Bclido and Kusenkos's model, the role of resonant field is 
played by the electromagnetic field interacting with this 
oscillatory field. 

We shall show below that the proposed model leads to 
a very efficient, although not thermal, photoproduction. 
We expect that in a second phase (quite the same as in 
Inflationary Cosmology) the production will thermalize 
with the surrounding matter and consequently changing 
the equation of state (EoS) of matter inside the star. Of 
course, a careful study of energy transport is needed in 
order to evaluate the impact of this kind of photopro- 
duction to the star evolution. In this work we just show 
that Parametric Resonance is a plausible mechanism for 
the pion-photon model leading to an anomalous and fast 
energy production in the interior of compact stars. The 
new, nevertheless important, questions above mentioned 
deserve, of course, a careful further study. 



II. PION CONDENSATION IN NEUTRON 
STARS 

Neutron stars have provided us with information on 
the properties of high density nuclear matter. Accu- 
mulated data have given useful information about reli- 
able EoS and thermal properties of the NS nuclear mat- 
ter which seem to suggest the appearing of some phase 
transitions there, e.g, pion condensation. The EoS then 
can be used to determine NS parameters, most notably 
the mass-radius relationship and the maximum allowable 
mass '2Q}. 

An interesting phase transition that can take place in- 
side a NS is the pion condensation. This possible conden- 
sation is due, to a large extent, to the high density of the 
neutron star medium. The possibility of pion condensa- 
tion in nuclear or neutron matter was initially proposed 
by Migdal [Tj , Sawyer and Scalapino Qj| . Whether 
pion condensate exist or not is of great interest since it 
can significantly enhance the cooling rates of NS [T^.ll5| . 
A very fast cooling for NS may be explained by a possi- 
ble pion condensation ^(| . PR mechanism start to work 
just after the condensate is formed as we shall show in 
the following sections. 

Condensed phase of negative pion (ir~ ) can take place 
if we neglect the strong pion correlations with ordinary 
matter in modifying the pion self energy. It is possible to 
show (see ref. [13) that the favorable condition to this 
assumption is satisfied for a neutron on the top of the 



Fermi sea to turn a proton and 7r when 

f^n l-^p fte > fTl-n 

where m w ~140 MeV is the tt~ rest mass. When this 
occurs, we say that the negative pion macroscopically 
occupy a single mode in its ground state, that is, they 
are in a condensate state. 

The 7T° condensation can take place also, since its ef- 
fective mass in the medium is zero 0]. Thus, the decay 
reaction occurs 



n— m + 7r . 

Possible condensation of positive pion (ir + ) will not 
taken into account here since Sawyer and Yao 01 had 
showed that in NS medium the number of ir + particles 
is much smaller than the ir~ particles. Therefore, even 
when 7r + condensation occurs it could be neglected when 
compared with the tt~ condensation. Neutral pion does 
not interact with the electromagnetic field of the NS and 
the 7r + condensate will does not affect the electromag- 
netic field in significant way, considering the above argu- 
ment. So, for the photoproduction, in our case, we can 
take into account only the interaction between the nega- 
tive pion condensate and the electromagnetic field of the 
NS. Henceforth, we will indistinctly use the wording pion 
condensate as 7r _ condensate. 



III. GENERAL FORMALISM IN THE SIGMA 
MODEL 

In order to get a tt~ condensate we will use the La- 
grangian formalism for the linear a model. Although the 
linear a model is a toy model it is enough, for our pur- 
poses, to do a qualitative study of the PR phenomenon. 

We start by considering the Lagrangian density for the 
electromagnetic field A^, the isoscalar a, the isotriplet 7?, 
and the nucleon isodoublet \t [l^. |20| : 

C = La + Ca + C 7r - 7T + + C n o + £y — Lint, (1) 

where 

C A = -^F^F^ + e 2 (A^ 2 , 
is the electromagnetic field contribution, 



C a = - (df.a) 2 + -m 2 a 2 - — cr 4 + cict, 

is the a contribution. Note that c\o term is responsible 
for the spontaneous symmetry breaking in this model. 
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represents the neutral pion contribution, 



C 



= d n TT <9 M 7T 



+ im 2 (2 7 r- 7 r+) - ^{2^i,+ f : 
is the charged pion contribution, 



1 



4r° = ^{d^f + 



denote the nucleon contribution and finally, the interac- 
tion Lagrangian density reads 



Ant = -ie(7T + 9 M 7T -7T 9 M 7T+)A^ - A" - (1 + T 3 ) + + iT • 7T75)]* 



A 1 
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[a 2 7r° 2 + a 2 (27r-7r + ) + 7r° 2 (27r-7r+)], 



r 



where Tj = 1,2,3 are the Pauli matrices and 



75 = 



1 

1 



The nucleon isodoublet '5 contains four-component 
proton and neutron spinor fields and can be written as 



and 7r denotes the pion field operator that annihilates a 
7r~ (or creates a which we can write as 



7Ti ± Z7T2 



(2) 



where 7Ti and W2 are real components of the pion field 
and if — (tt , 7r + , 7r~ ) . 

The constants in the Lagrangian density can be related 
(in tree approximation) to the physical quantities /7r~93 
MeV, MeV jHlH: 



2A 2 / 7r 2 = TO CT 2 -TO 7r 2 , 



2m 



Note that if we add a symmetry breaking term in the tt 
Lagrangian density, instead in a Lagrangian, then differ- 
ent constants must be defined |2l| . Finally e = is the 
electromagnetic coupling constant and g — J^rra, where 
in is the nucleon mass (typically m~940 MeV). Usu- 
ally the mass of a meson is taken as bigger than 1 GeV 
[Mil]. Without loss of generality, we take m a = 1020 
MeV, and then we get to -700 MeV, A~7.7 and #~87420 



MeV . A priori we do not have information about the 
a meson and in the NS medium we have a lot of candi- 
dates to be the "er meson" [l(| • Since this is a qualitative 
study we choose this value for m CT due the mass of <f> me- 
son, whose mass is around 1019 MeV. 

Harrington and Shepard 0, |2(j, using (QJ, showed 
that the pion condensate and the strong electromagnetic 
field of a NS may co-exist in a particular region (shell) in- 
side of the compact star, as well nucleons. They showed 
also that the pion condensate can assume a supercon- 
ducting behavior, even in the presence of an intense elec- 
tromagnetic field. According, we assume in this work 
that the pion condensate and the electromagnetic field 
may co-exist inside of the NS, as showed by the authors 
of 0) l20j . In terms of Harrington and Shepard's work 
the pion condensation, without a superconducting phase 
transition, occurs when the neutron star magnetic field 
H is weaker than some critical low value, H c \. This sets 
the basis for our physical model. 

Now, Parametric Resonance phenomenon belongs to 
a class of ordinary linear second order differential equa- 
tion called Hill's Equation This class possess two 
distinct important cases (but they are not the only ones, 
of course): the Mathieu's and Lame's ones. In the con- 
text of particle production by the PR mechanism, Math- 
ieu's Equation was firstly used by Traschen and Branden- 
berger y| and Dolgov and Kirilova Q and Lame's Equa- 
tion by Kofman et al Q . Following the above authors we 
find the equation of motion (EoM) for the pion conden- 
sate and its periodic solutions which, in turn, enter as 
periodic potentials in the EoM for the electromagnetic 
field. These are Hill's Equations whose the most use- 
ful ones are the Mathieu and Lame types as mentioned 
above. 

The specific EoM for A^ depends only on the form 
of periodic solution for the pion condensate EoM. If its 
periodic solution is a simple trigonometric function then 
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the A^'s EoM is classified as Mathieu's type; if its so- 
lution is an elliptic function, the EoM is classified as 
Lame's type. For both types, solutions for A M 's EoM are 
guaranteed by a Floquet's theorem [(| and the behavior 
(growth) of a particular solution is determined, basically, 
by its Floquet exponent. If the Floquet exponent is a 
real number, the solutions are unstable, i.e., they present 



resonant bands in their parameter space in which the so- 
lutions growth (or decay) exponentially. In the case of 
imaginary number, the solutions are stable (bounded). 

So, in order to find the EoM for the electromagnetic 
field we must find firstly the EoM to the tt~ condensate, 
applying the Euler-Lagrange equations in (JJJ. We get, 



□tT 



-m 



+ A 2 7t + (t0 2 + i 9 V2*[r_7 5 ]1' = 0, 



(3) 



r 



where t_ = (n — «T2)/\/2 and the back reaction of 
field over 7r~ condensate was neglected. As mentioned 
above, our purpose is to study the effects of 7r~ con- 
densate over the A^ field and then investigate the PR 
phenomenon inside of an isolated NS. Back reaction of 
A^ field can affect the evolution of ir~ condensate, but 
we expect, in a first order approximation, that this con- 
tribution may be neglected as in the preheating phase of 
the Inflationary Cosmology 0,0 . 

The 7T° field does not interact with the A^ field, so we 
can take it in its ground state 







tion ©. 

Symmetry breaking term presents in C a (equation Q ) 
leads to a a condensation prevent us to consider (o) vac = 
pll |22| . In order to estimate the a contribution we 
take (a) = b(a)y f , r ., where b is a positive real number. 
Following [l9ll20ll2lll2^| we may write 



(<r)vac—U cos 6 

where 0^0 is the chiral symmetry angle that gives infor- 
mation about the a and 7r condensation. So, 



and this implies a non condensation phase to ir° meson 
The 7T~ condensate is not affected by this assump- 
tion since conditions on both condensations are inde- 
pendent, as seen above. Now taking the approximation 
(7r° 2 )~(7r ) 2 we can neglect the tt° contribution in equa- 



= b( 



-afn, 



where a = bcos9 > is a real parameter. 

Then, equation Q reads, after some algebra with the 
nucleon components, 



J 



□tT 



2 „2 



m 



n~ + A 2 7r+(7r-) 2 
I 



>/2|*r. 



0. 



(4) 



If we take values previously obtained for too, A and /„•, 
we obtains two distinct cases for equation J1J: 



case +2M 2 to a 2 > 1.92. So we can write equation (@J 
as 



a 2 > 2m ° 2 



d.92, 
d.92. 



For the sake of simplicity we can take into account 
these two inequalities defining 



\' 2 f 2 n 2 

A U - too 2 = ±2M 2 , 



where we have case — 2M 2 corresponds to a 2 < 1.92 and 



□tt" ± 2Af 2 7r~ + A 2 7r + (vr") 2 - igV2\^ n \ 2 = 0. (5) 

Also, as a first approximation, we can consider here 
a region inside the NS in which 7r~ condensate has a 
homogeneous and isotropic distribution. So, we can ne- 
glect spatial contributions and then equation (J3J) may be 
rewritten as 



dt 2 



7T~ ±2M 2 7T- + A 2 7T+(70 



0. (6) 
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We will analyze the above equation (with ± signs) in 
two cases: Homogeneous (\^ n \ — 0) and Inhomogeneous 
(l^nl 7^0) equation. We will see that for both signs we 
get periodic solutions which lead to PR for the electro- 
magnetic field . 

A. Homogeneous Case ( |^n| 2 = 0) 

Neglecting nucleon contributions we can write © as 

-^U-±2M 2 7r-+A 2 7r+(O 2 =0. (?) 
at 2 - 

This ordinary nonlinear second order differential equa- 
tion is not in any way a simple one. ft can be solved 
in terms of Weierstrass's Elliptic Functions [25| which 
are periodic. We can use relation (J2J) and write 10 as a 
real ordinary nonlinear second order differential coupled 
system: 

^TTi ± 2M 2 7T! + ^TTi 3 + ^-^TT, = (g) 

±zir 2 ± 2M 2 ir 2 + ^ir 2 3 + ^-7ri 2 7r 2 = 0, 

where the same sign must be taken simultaneously in 
both equations for the factor 2M 2 . 

It is easy to see that system described by © is sym- 
metric under 7Ti< — m 2 exchange. This fact suggests sim- 
ilarities between amplitudes, wavelength and frequencies 
of 7r~ and 7r + solutions. In the next subsection we will 
make a numerical analysis of this system. 



1. Numerical Solution 

Figure shows the evolution of the system (jSJ) for the 
case -2M 2 with M = 140 MeV and M = 420 MeV. 
These values were obtained when a 2 = 1.78 and a 2 = 
0.54, respectively. Since a 2 < 1.92, these values can be 
completely arbitrary (recall, a > always). We took 
some values for M as one pion mass and three pion mass, 
approximately, just for numerical analysis. A priori, any 
real positive value for Af may be used. 

Figure [21 shows the temporal evolution of the system 
© for the case +2 A/ 2 for M = 140 MeV and M = 420 
MeV. These values for M were obtained when a 2 =2.1 
and a 2 = 3.3, respectively. 

In both Figures, n and |3 the auxiliary 7Ti and n 2 have 
a periodic behavior. Therefore, as a direct consequence 
of this fact, the pion condensate will have a periodic be- 
havior by 0. But it is still necessary to know what kind 
of regularity we have: a trigonometric one or an elliptic 
one. 

Figure |3 shows the phase space (7Ti,2, ^1,2) for the case 
— 2Af 2 . Figure 0] shows the phase space for (711^2,^1.2) 
for the case +2A/ 2 . Both pictures confirm the periodic 



Case 


TV IT Ti T t 7 

Mass MeV 


Wavelength (I) 


frequency (v) 


(-2M 2 ) 


M = 140 


~ 0.1499 


~ 6.6675 


(-2M 2 ) 


M = 420 


~ 0.0541 


~ 18.4942 


(+2M 2 ) 
(+2M 2 ) 


M = 140 


~ 0.0317 


~ 31.5111 


M = 420 


~ 0.0106 


~ 94.5337 



TABLE I: Table with wavelength and frequency for M = 140 
MeV and M = 420 MeV in -2M 2 and +2M 2 cases. 

behavior of solutions for tti and n 2 and suggest us that 
periodic solutions shown in Figure^and[5]have an elliptic 
function behavior, that is, they could be described by Ja- 
cobi's Elliptic Functions. Therefore, as we will see, Hill's 
Equation that comes out from for the ^4^'s EoM is, 
in fact, a Lame's Equation whose solutions are of elliptic 
type. 

It must be observed that JHJl is an ordinary second or- 
der differential nonlinear coupled system and this may 
imply a strong dependence on a given set of initial con- 
ditions. Alongside this paper we adopt, without loss of 
generality, 7Ti(0) = and 7Ti(0) = 1 (the same values for 
ir 2 , of course). 

It is possible to obtain, numerically, an estimative of 
wavelength (/) and frequency (v) for the cases — 2M 2 and 
+2Af 2 , if we consider roughly v^d~ x . The Table I show 
some values of I and v. 

Particle production by the PR mechanism depends 
strongly on the amplitude and frequency oscillation of the 
source field, that is, on the behavior of the ir~ field. De- 
termination of these parameters provide a roughly quan- 
titative estimate of the particles number produced in the 
NS medium. 



2. Exact Solutions 

The symmetry of system © and its numerical analysis 
suggest us the ansatz 

7T2 = TTl , 

and the system is now reduced to a single ordinary non- 
linear differential equation written as 

d 2 

— TTl ± 2Af 2 7Tl + A 2 tti 3 = 0. (9) 
dt 

The above equation is the so-called Duffing's Equa- 
tion [2^| and can be solved in terms of Jacobi's Elliptic 
Functions. This fact is in complete agreement with the 
above numerical analysis. The Duffing's theory allow us 
to understand the behavior of exact solutions as follow 
below. 

Equation 10 can be transformed in a following system: 
I $ =-W±2M 2 n 1 , 
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60 
40 
20 

"l,2 
-20 
-40 
-60 



0.1 0.2 0.3 0.4 0.5 

t 



FIG. 1: For the case -2M 2 : (a) M = 140 MeV and (b) M = 420 MeV. Note that 7Ti and ty 2 are similar in (a) and (b). In 
both cases, A = 7.7. 
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(b) 



. 004 

0.002 
"1,2 0- 
-0.002 
-0 . 004 



).l 0.15 

t 



. 0015 
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-0.001 
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FIG. 2: For the case +2M 2 : (a) M = 140 MeV and (b) M = 420 MeV. In both cases, A~7.7. Again, they are similar in (a) 
and (b). 



where rj = r/(t) is a real function. The field direction of 
this system in the phase plane (77, 7Ti ) can be written as 

dtj dt _ drj A 2 vr 1 3 ±2M 2 7r 1 
dt diTi dixx rj 

Direct integration of this equation results 
A 2 

V 2 + y ""I 4 ± 2M 2 tti 2 = 2c, (10) 

where c is an integration constant with dimensions of 
MeV 4 . Figure|S]show a pictorial representation of {rj, iti) 
phase space. In Duffing's theory context, this c parame- 
ter can be interpreted as an energy level to a given closed 
trajectory. For the case rj 2 + +2M 2 iti 2 = 2c, we get 
the so-called hard spring case. The solutions are closed 
trajectories centered at the origin of the coordinate sys- 
tem, as showed in Figure 01. In this case, only positive 
values for c are permitted. For rj 2 + ^iri 4 — 2M 2 7Ti 2 = 2c, 
we have the soft spring case and, like the previous one, 
the closed trajectories, which exist only near the center of 
the potential, represents yet periodic solutions. For these 



situations, each value of c represents an energy level for 
a given trajectory, as showed in Figure \5]p. In this case, 
far from the center of the potential, the solutions are 
non-periodic, for an arbitrary but fixed value of c |26|. 

For the hard spring case we have only one real min- 
imum at (77 = 0, 7Ti = 0). For the soft spring case the 
center of the coordinate system is an instability point 
given by {rj = ±^/c,tti = 0) which separates two max- 
ima, each one localized near of the system center and are 
given by m = ±&M. [2^. 

A separation of variables followed by an integration 
take equation (JTOJ , for both cases (+2M 2 and -2M 2 ), 
into the form: 



r z dv 1 - 

\ ; = = V2Mt, (11) 

where, for the sake of simplicity, we took t = 0, and 
defined z — -^^\ and f = 

The above integral depends strongly on the values of 
p. So we split it in several cases as showed below. 
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FIG. 3: Phase space for the case -2M 2 : (a) M = 140 MeV and (b) M = 420 MeV. In both cases, A = 7.7. The shape of these 
curve suggests periodic, bounded solutions as elliptic functions, as can be seen in [2^|. 





015-0.001-0.0005 0.0005 0.001 0.0015 



FIG. 4: Phase space for the case +2M 2 : (a) M = 140 MeV and (b) M = 420 MeV. In both cases, A = 7.7. Also the elliptic 
behavior is suggested by these graphs. 



Exact Solutions: Case —2M 



p < — 1/2: In this case we do not have real solutions 
for the integral presents in iJUJ. 

p > 0: It is possible to write the roots of — z' 4 — z' 2 + 
f = as 



1 1 



1 1 



r 2 = + 2 V 71 + s 2 = +- + -y/l + 2p, 

and using the table of integrals in [23 we write the inte- 
gral in fTT)l as 



I z r 2 + s 2 \ o 



and q 2 is the elliptic function modulus (0 ^ q ^ 1). 
From (fill) and i|12[l and after some algebra one obtains 



A 



2psn 



2M" A v / x /TT2^(2 + sn 2 )-sn 2 ' 
where, for short, 



(13) 



dz' 



-;Kh,q% (12) 



o v/(z' 2 + r 2 )(s 2 -z' 2 ) Vf^T? 



where K[y,q 2 ] is the elliptic integral of first kind given 
by 



ch, 



° v /(l- (Z 2 )(l-sin 2 7 )' 



and 



[4(l + 2p)]«Aft, 



1 + VTT2p 



is the Jacobi's elliptic sine that possess a double period. 
The behavior of this solution can be fixed by an appro- 
priate choice of the parameters M and p. 

— l/2<p<0: Again, using [23] we obtain two distinct 
solutions which depends on the roots 



1 1 



1 1 



+ 77+ oVl + 2p, * 2 = +---Vl + 2p, 
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7Ti 




7Tl 




FIG. 5: Pictorial representation of (a) hard spring case and (b) soft spring case. Note that in the soft spring case we have two 
instability points 7ri = ± V ^ M . 



Using [23 we can write 



Exact Solutions: Case +2M 



dz' 



y/-z'*±z'* + f 



where now 



-K[C,q% 



(14) 



In the case +2M 2 , using [27]], we obtain just only one 
solution. 



p > 0: We can write the roots of — z' 4 + z' 2 + | — 



as 



£ = arcsin 



r / z — s £ 



Again, from equations ((14(1 and 1)11(1 we get 



A = 

2M 7 " 1 A /l + v /l + 2p(l-2sn*)' 



(15) 



Using [27J and writing the roots as done in the case 
above one obtains 



^/-z /4 ±z 



dz 1 , 
= = -K\ip,q 1 , 

2 



(16) 



where 



1 1 



1 1 



and we obtain the same solution expressed by equation 
(fH!j) . Solutions obtained until now does not possess sin- 
gularities and have a time periodic evolution. Jacobian 
sine behaves, roughly speaking, as a trigonometric sine. 

In our numerical analysis we noted that for an appro- 
priate fixed M (140 MeV or 420 MeV), frequencies os- 
cillation versus p of equations in ((13(1 , (|15fl , (|17fl presents 
a logarithm form, that is, for some fixed value of p, the 
frequency have a slow, but growth behavior. The PR 
phenomenon here depends on frequency oscillation of tt~ 
field and it indicates that, from a determined value of p, 
the particle creation rate itself have a slow, but continu- 
ous growth. 



ip = arcsm 



1 2 = 



Finally, from equations (fT?)|l and (fTT]) one obtains 



2M* 1 = V 1 + V /l + 2p(l-2 S n 2 ), (17) 



where, in both cases, 



2VTT2^ 



1 + Jl + 2pMt, — v , F 
V y ' 1 + ^TT2p 



Again, all characteristics of these solutions are con- 
trolled only by the free parameters M and p. 



B. Inhomogeneous Case 7^0) 

The l^&ril 2 field will be treated here, for the sake of 
simplicity, as an isotropic, homogeneous medium where 
interactions between the pion condensate and the A 11 field 
can take place. In other words, |\& n | 2 — ^n^n possess a 
constant value in the NS medium. 

Here, \^> n \ 2 describe a neutron matter in its condensate 
phase. The condensation of neutrons can take place in 
the usual way, the neutrons comprise Cooper pairs which 
obey the Bose-Einstein statistics. In the condensed phase 
its ground state assume a non-zero value, that is 



^0. 
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In this phase it is well-known that neutrons can be 
described by a standing wave function. We take the ap- 
proximation fy n = ($ n ) and write 

where k is the momentum of neutron condensate. If we 
assume isotropy of the condensate neutron we can write 



where g is the coupling constant presents in the interac- 
tion Lagrangian. In order to get a qualitative study we 
take h = g(^n) vac as a free parameter and analyze the 
behavior of the pion condensate for some values of h. 

Virtually, the above system does not have an exact so- 
lution. So, we can make only a numerical study of its 
temporal evolution. Figure shows the temporal behav- 
ior of 7Ti (a) and tt 2 (b) for (— 2M 2 ). Note that regularity 
persists even for a wide range of h; Figure shows the 
time behavior of 7Ti (a) and TT2 (b) for (+2M 2 ). Again, 
the regularity persists. Some irregularities in the figures 
are due to numerical approximation. 

In this paper we will only analyze below the PR phe- 
nomenon through exact solutions of the pion condensate 



Note that the interaction term 2e 2 (7r~7r + )v4/ i is the 
most interesting in the above equation to our purposes. 
As we will see in a moment below, it is responsible for 
the PR phenomenon. We need to make some approxima- 
tions in order to get an EoM which can be handled more 
easily. We start by neglecting nucleon contributions, i.e., 
we firstly study the Homogeneous Case ( hfr J 2 = 0). This 
can be made to first order approach |19l l2Cj . In the next 
subsection we study this case. 



A. Photoproduction in the Homogeneous Case 

(l*n| 2 =0) 



it as 



(*„) - {Vn) vac exp[ikr], 

Taking into account this result and using the relation 
(|2"|) we obtain a nonlinear second order non-homogeneous 
system 



= -.g(*, i ) 2 , QC sin(2fcr) 
= -3<* n ) 2 t , QC cos(2fcr), 

I 

EoM. The fact that inhomogeneous system (|18|) possess 
periodic solutions suggests that the PR phenomenon can 
persists even in this situation. 



IV. PHOTOPRODUCTION BY THE A" FIELD 
OSCILLATIONS 

In this section we study solutions obtained for the time 
dependent electromagnetic field coupled to the pion con- 
densate field. From the Lagrangian density (JTJ one ob- 
tains the EoM of the electromagnetic field which reads: 



(19) 

I 



OA" - d»d v A v + ie(w + d fJ "K- - n~d^n + ) + 

+2e 2 (7r~ K + )A f * = 0. 

Now, from the symmetry of pion system in (8) we can 
assume, for the sake of simplicity, the ansatz -ki = m the 
term ie(ir + 'd^ir - — n~d ,J "7r + ) vanishes. This ansatz was 
suggested by the analysis done above of the solutions for 
the pion condensate solutions obtained from the homo- 
geneous system for the case l^nl 2 = 0. Thus, the above 
equation reduces to 

DA^ 1 -d^d u A u + 2e 2 (n-TT + )A^ = 0. (20) 
Introducing the gauge field transformation 



I^ttx ± 2M 2 tt 1 + l-TTj 3 + ^-7T 2 2 7r 1 

^TT 2 ± 2M 2 1T 2 + \^ + 



I 

UA^ - d»d v A v + iein+d^n- - n-d>*n + ) + 2e 2 {n- tt+)A» - e^-V(l + r 3 )* = 0. 



For l^nl 2 = the equation i|19|) reads 



(21) 
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FIG. 6: Figures above show the behavior of ix\ and TV2 for the case —2M 2 . For (a) and (b) M = 140 MeV: continuous- line 
(h = 1000) and dashed-line (h = 0.001). For (c) and (d) M = 420 MeV: continuous- line (h = 1000) and dashed-line (h = 0.001). 
Note that even under extreme variation of h the regularity still persists (h — g{^n va c} 2 )- 




FIG. 7: Figures above show the behavior of m and 1V2 for the case +2M 2 . For (a) and (b) M — 140 MeV: continuous- line 
(h = 1000) and dashed-line (h = 0.001). For (c) and (d) M = 420 MeV: continuous- line (h = 1000) and dashed-line (h = 0.001). 
Note that even under extreme variation of h the regularity still persists (h — g(^n V ac) 2 )- 



where £(r, t) is an arbitrary function of r and t and using 
the Lorentz gauge condition 



since —d^d v A v vanishes. An isotropic gauge field A"* 
can be written as 



A'^x) = X (x)e% (22) 

where xi x ) can be expanded in Fourier modes and is 
a polarization vector. 

Substituting (|22|l in l|20|) and taking into account the 
Lorentz gauge condition one obtains 
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e"px(as) + 2e 2 (7r-7r+)x(a;)] = 0. (23) 

Now, in a crude, but nevertheless unrealistic approx- 
imation, we can take also the electromagnetic field as 
homogeneous inside the NS. So l|23[l reduces to an ordi- 
nary differential equation in the time variable. Taking 
the Fourier transform to k space momentum we finally 
get 

^Xk(t) + [k 2 + 2e 2 (n-n+)] Xk (t) = 0. (24) 

This is a homogeneous differential equation for an os- 
cillator with frequency (or mass) varying in time. We 
define the frequency oscillation of the Xk(t) ne ld as 

uj k = k 2 + 2e 2 (7r _ 7r + ). 

For the Homogeneous Case (|\? n | — 0) and solutions 
of 7r~ and 7r + fields, previously obtained, we may write 
{^fy as a Lame's type equation. 

A general Lame's Equation have the form 

^i/>(t) + [A + BQ(t)Mt) = 0, 

where A and B are real numbers and Q(t) is an elliptic 
function. Phase space of the parameters A and B is not 
easily defined, and this is exactly the case of equation 
(|24p. Nevertheless it still possess an important property 
namely, the existence of instability (exponential growth) 

Xfe(*)~exp[/ifct], (25) 

where [i k G (— oo,oo) is the so-called Floquet exponent 
(or characteristic exponent) labelled by the k momentum. 
To find exact solutions for the Lame's Equation is a very 



where 

_ 2 16M 2 e 2 psn 2 
UJk= + A 2 [VTT2^(2 + .sn 2 )-sn 2 ]' 



difficult task. For example, if the coefficient of (7r~7r + ) 
term in (|24[1 may be written as n{n 4- l)q 2 , where n is an 
integer and O^q^l is the elliptic function modulus, then 
an exact solution can be founded. On the other hand for 
a qualitative study we can use [i k taken from the formal 
solution of Hill's Equation [f| , which can be written as 

Mfe = I cosh -i[l_A( ) sm 2 (^)], (26) 

7T 2 

where A(0) is an infinite determinant. Figure [S] shows 
[ik > for arbitrary values of k > momentum with 
several fixed A(0) values. 

The Floquet 's theorem @ gives 

Xk(t) = exp[[i k t]p(t), 
where pit) = p(t + T) is a real limited periodic function 
with period T. When [i k > 0, instability corresponds to 
an exponential growth of occupation number 

n k (t)~exp[2n k t], (27) 

that may be interpreted as number density of particles 
produced by PR mechanism 0, Q| • 

In the following subsections will analyze the equation 
by means solutions ifP3)l. (|T3|) and (JT^ for the ir~ 
condensate. 

1. Homogeneous (\^„\ 2 = 0): Case -2M 2 

In order to obtain the Lame's Equations, the solutions 
obtained for the pion condensate EoM were substituted 
in equation l|24l) . All solutions that we found for the non- 
linear field (source) are periodic and substituting them in 
(12411 we got a Hill's Equation. For the sake of simplicity, 
below we show some of them. 

Using equation l|13|) and after some simplifications we 
write P3)l as 



Xk(t) = 0, (28) 



I 

having a Jacobian periodicity, and, a priori, may be write 
in a more handle way, i.e., as just one Jacobian sine or a 
sum of them. For our purposes, it is sufficient in the Hill's 
theory that the potential (|13|l has a periodic behavior. 



alt 2 



Xk(t) 



16MVp 



A 2 [yr+2^(2- 



■rSfl 



The above equation may be classified as a Lame's type 



Using solution (|15|) . equation (|24|l can be written as 
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FIG. 8: Floquet exponent fit as a function of k and fixed values of the determinant A(0): (a) A(0) = —1000, (b) A(0) = —100, 
(c) A(0) = -1 and (d) A(0) = -0.001. If A(0) < 0, fi k is real; if A(0) > 0, fi k is imaginary. 



dt 2 



Xk{t) 



16MVp 



A 2 [l + x /l + 2p(l-2sn 2 )] 
I 



Xk(t)=0, 



(29) 



where clearly 



ojk = k + 



16MVp 



A 2 [l + vT+2p(l - 2sn 2 )]' 



We got again an equation whose potential is written in 
terms of elliptic functions and it is of Lame's type. 

Now, using (fT7|) . equation (|2"l|) can be written as 



dt 2 



Xk(t) 



, 2 8M 2 e 2 (l + ^/T+2p) 
C + X 2 



16M 2 eVl + 2 P 
A 2 ' 



Xk(t) = 0, 



(30) 



with 



2. Homogeneous (\^ n \ = 0): Case +2M 2 



j2 8M 2 e 2 (l + y'l + 'Ip) 16A/ 2 e 2 Vl + 2p 2 
u>k = k H — — sn 



A 2 



A 2 



This is, clearly, a Lame's Equation which is easier to 
handle than the another ones above. 



In this case we get just one equation. Using $Tj^ equa- 
tion (|24(l may be written as 



J 



d 2 , . 
^Xk(t) 



16M 2 e 2 p 



X 2 yi + 2p(2 + sn 2 ) - sn 2 } 
I 



Xk(t) = 0, 



(31) 



where 



16MVp 



X 2 yi + 2p(2 + sn 2 ) - sn 2 
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represents its mass oscillation coefficient. 

All Lame's Equations obtained above have only two 
free parameters, M and p. These two free parameters 
control frequency and amplitude of pion oscillations and 
consequently by equation (|27|l also the rate of particle 
production of the photon field. 

Energy density for each Lame's Equation can be cal- 
culated by 

1 f kc 

p 1 = — J dkL0 k n k (t)k 2 . (32) 

We are not taking into account here the back reaction 
of the particles produced. 



Now we take as zeroth approximation the nucleons con- 
tribution as a constant in space and time variables. This 
is, of course, a crude approximation relied on isotropy, 
homogeneity and independence on time for the nucleon 
matter inside the star core. Nevertheless the figures El 
and [7| shows that periodicity of the pion fields is robust 
against extreme variations of the nucleons contributions 
(in zeroth order). That is, the nucleons have only a weak 
influence on the pion evolution, whereas the strong influ- 
ence on amplitude and frequency comes from the value 
and signal of pion mass term M. In adition we believe 
that so restrict conditions of Constance of nucleons con- 
tribution can be relaxed and the periodicity of pion fields 
will not strongly affected. This is enough in order to PR 
mechanism works well. Therefore, for our qualitative ap- 
proach we can neglect the nucleon contribution. Taking 
the Fourier Transform, for the approximation described 
above, we obtain the same Lame Equation (24). 

All numerical solutions suggest that equation (24) pos- 
sess resonant bands, that is, the PR phenomenon may 
occurs even in the presence of nucleons. Clearly, since n k 
has exponential behavior, so p k in 1)32(1. 



V. PARAMETRIC RESONANCE AND 
NEUTRONS STARS FLARES 

We had showed that Parametric Resonance phe- 
nomenon, under certain conditions, may occur inside of a 
NS modifying briefly its thermal evolution. It is worth to 
mention that our approach can be implemented, not only 
for the Harrington and Shepard's model, but in fact for 
any NS model which presents superconductor phase with 
periodic motions (EoM) and a suitable coupling to Elec- 
tromagnetic Field to furnish a Hill's Equation. PR mech- 



In order to obtain a more realistic physical descrip- 
tion is necessary introduce by hand a term, in Q or, 
yet in l(24|l , that represents the back reaction of particles 
produced. Since this work aims just to propose a new 
process of photoproduction in NS, it is sufficient to prove 
that the PR phenomenon could exist in the NS medium, 
as it is shown above. 



B. Photoproduction in the Inhomogeneous Case 

(l*n|V0) 

Using the gauge transformation (|21|) with the Lorentz 
condition, equation (|19|) can be written as 



(33) 

I 

anism can be effective enough to produce a huge num- 
ber of photons, out of equilibrium with the surrounding 
matter, which must diffuse out from the NS. In this case, 
clearly, a subsequent reheating of the NS interior must 
occur, much the same as in the Inflationary Theory, but 
in a far smaller energy scale. From the mathemartical 
point of view we have two cases to consider, namely, for 
real Floquet exponent and when it is an imaginary num- 
ber. Both cases lead to an enhancing of the NS bright- 
ness, but the possible observable effects are expected to 
be very different as detailed below. 

Taking for granted the current theories of NS interi- 
ors as nearly correct, in our opinion, this kind of PR 
phenomenon can, in principle, be detected by observa- 
tion of anomalous enhancing of the NS brightness. The 
intensity, as well the lapse of time of this higher bright- 
ness depends on the details of diffusion of the produced 
photons to reach the surface of the star. This is a diffi- 
cult problem since, by its turn, it depends on the model 
of the NS matter structure. On the other hand, recent 
works seem to indicate a nearly black body surface ra- 
diation from Neutron Stars whose details also depends 
on the model of their atmospheres [30l l3lL I32I l33 |. If 
this is correct, in the case that PR mechanism is active, 
it leads to some interesting results on the radiation re- 
leased by a NS. We consider first the case of real Floquet 
exponent. On a blackbody spectrum background, PR 
manifests its presence, by the theory outlined above, in 
a weak regime, as a perturbation on NS spectrum with 
superimposed bumps on the equally spaced frequencies 
due to the periodic pattern of the Floquet exponent as a 
function of the momentum k, as is clearly seen in equa- 
tion i|26|) and in Fig. [S] This periodicity implies that high 
energy photons are also produced and their spectral lines 
should dominate over the blackbody ones, which lead us 



I 

UM 1 + iein+d^TT- - 7r"a M 7r + ) + 2e 2 (tt - tt + ) - e^~$ (1 + t 3 )# = 0. 



14 



to expect that PR effects could be better measured on the 
tails of a blackbody spectrum. Furthermore, although in 
a qualitative ground, in the case of a very strong regime 
(with a large real Floquet exponent) of photon produc- 
tion, the external star layers can not resist against the 
enormous pressure due to the reheating of matter in the 
star core by high energy photons. In this case the star 
ejects its upper layers with a huge burst of energy releas- 
ing. In this way our model suggest a feasible and simple 
explanation for the recent observations of neutron stars 
giant flares 0, l35l l3(| . Shortly, PR mechanism with a 
real and large Floquet exponent is better to explain non- 
periodic, but giant flares from NS. On the other hand, if 
the Floquet exponent is imaginary, Eq. (27) shows that 
the number density of produced particles is a periodic 
function in time and so we can expect a periodic pertu- 
bation on the brightness of the NS while the PR regime is 
active. This is quite interesting since giant periodic flares 
from repeaters have been observed in the past few years. 
See, for example, |36|. However, for this case, we can 
have only an incomplete view since for giant flares the 
unkown coefficient in Eq.(27) must be very large, since 
the time factor is a trigonometrical function and so a 
bounded function. So, a further study is needed to solve 
this quest. 

VI. CONCLUSION 

Firstly, it must be stressed that Parametric Resonance 
in NS can be a rare phenomenon, due to the special condi- 



tions of the superconductor matter and its coupling with 
the electromagnetic field. Nevertheless, in the weak PR 
regime, as mentioned above, it is worth to notice that, 
in principle, it is possible to measure an uprising regular 
(equally spaced frequencies) brightness variation on the 
nearly blackbody spectra, perhaps better measured on 
their tails, of those neutron stars close enough to us. In 
adition, in the strong PR regime, the model offers also 
a simple explanation for some reported rare giant flares 
coming from these compact objects. On the other hand 
the signal can also presents a time periodic pattern. In 
this case, as in Inflationary Theory, it is not expected an 
explosive behaviour of the star, but only periodic flares 
with low intensity. 

Finally, more importantly perhaps, we are inclined to 
think that Parametric Resonance is a phenomenon active 
not only in Inflationary Cosmology but in several much 
smaller coupled physical systems with extreme conditions 
of pressure and temperature in our Universe. 
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